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Abstract 

We study the baryonic charmonium decays of B mesons, B + — ► rj c K + and B + — > J/ipK + , where 
7] c and J/ip subsequently decay into a pp or AA pair. The charmonium produced in the above B 
meson decays is fully polarized. The polar angular distributions of the baryon-antibaryon pairs 
are presented, along with fit results to a 1 + ob cos 2 # parametrization. Comparisons are made 
with the results from e + e~ — > J ftp formation experiments. We also report the first observation of 
rj c -> AA. The measured branching fraction is B(rj c -> AA) = (0.871^ ±0.14 ±0.27) x 10 -3 . This 
study is based on a 357 fb^ 1 data sample recorded on the T(4S) resonance with the Belle detector 
at the KEKB asymmetric-energy e + e~ collider. 

PACS numbers: 13.25.Gv, 14.40. Gx, 13.40.Hq 
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There have been many reported observations of baryonic three-body B decays in recent 
years jH 0, 0, 0]. An interesting feature of these observations is the presence of a peak 
near threshold in the mass spectra of the baryon-antibaryon pair. Studies show that these 
enhancements are not likely to be resonance states, as the baryon angular distributions 
are not symmetric in their respective helicity frames 0. Other visible structures in the 
mass spectra arise from charmonium decays. It is natural to compare the baryon angular 
distributions from charmonium decays with those in the region of the threshold enhancement. 
There is a particular interest in J ftp — > pp, where the proton an gula r distribution has been 
studied by DASP @, DM2 @, MarkI 0, Markll and BES jD). W mesons produced 



in e + e~ —>■ J/ip are transversely polarized. Accordingly, the baryon angular distribution can 
be parameterized as ~ 1 + a cos 2 6, where 9 is the baryon polar angle in the J/ip helicity 
frame. Many theoretical predictions [l2| exist for the value of a. The current world average 
value of a, obtained with above measuremens from J/ip —> pp decays, is 0.66 ± 0.05. 

The study of two-body baryonic decays of charmonia at a I?-factory has several different 
features comparing with an e + e~ machine running at the J/tj) mass. J ftp mesons from the 
decay of spinless B mesons are fully longitudinally polarized. This provides a useful cross 
check for previous measurements with transversely polarized J/ip mesons. The charmonia 
from B decays do not suffer from the beam hole effect, such that events with | cos#| near 
1 can be detected. These events are effective to determine a. A 5-factory is also immune 
from e + e~ — > qq — > pp background, where q stands for auorrf quark. For previous studies 
this background is intrinsically embedded and hard to separate on an event-by-event basis. 

We use a 357 fb _1 data sample, consisting of 386 x 10 6 BB pairs, collected by the Belle 
detector at the KEKB asymmetric energy e + e _ (3.5 on 8 GeV) collider [13]. The Belle 
detector is a large solid angle magnetic spectrometer that consists of a three layer silicon 
vertex detector (SVD), a 50 layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time of flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of Csl (Tl) crystals located inside 
a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux return 
located outside of the coil is instrumented to detect K® mesons and to identify muons. The 
detector is described in detail elsewhere 3- 



In this study of two-body baryonic decays of charmonia we focus on the decay processes, 
B + — > ppK + and B + — ► KKK + ([15:]. The event selection criteria are based on informa- 
tion obtained from the tracking system (SVD+CDC) and the hadron identification system 
(CDC+ACC+TOF). All primary charged tracks are required to satisfy track quality criteria 
based on the track impact parameters relative to the interaction point (IP). The deviations 
from the IP position are required to be within ±1 cm in the transverse (x-y) plane, and 
within ±3 cm in the z direction, where the z axis is defined by the positron beam line. 
Proton, kaon and pion candidates are selected using p/K/n likelihood functions obtained 
from the hadron identification system. For the primary protons from B decays we require 
L p /(L p + L K ) > 0.6 and L p /(L p + L n ) > 0.6, where Lp/x/^ stands for the proton/kaon/pion 
likelihood. We require Lk/{Lk + L n ) > 0.3 to identify kaons. A candidates are recon- 
structed from decays into the pir~ channel. Each candidate must have a displaced vertex 
and flight direction consistent with a A originating from the interaction point [16]. To reduce 
background, a L p / (L p + L n ) > 0.6 requirement is applied to the secondary proton from the 
A decay. 

To identify the reconstructed B meson candidates we use the beam energy constrained 



mass, M bc = J E£ cam - p 2 B , and the energy difference, AE = E B - E heam , where E heam is 
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the beam energy, and ps and Eb are the momentum and energy of the reconstructed B 
meson in the rest frame of the T(4S). The candidate region is defined as 5.2 GeV/c 2 < 
M bc < 5.29 GeV/c 2 and -0.1 GeV < AE < 0.2 GeV. From a GEANT based Monte Carlo 
(MC) simulation, the signal is peaked in the region 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 and 
\AE\ < 0.05 GeV. 

The dominant background arises from the continuum e + e~ — > qq process. The back- 
ground from b — > c and charmless mesonic decays is negligible. In the T(4S) rest frame, 
continuum events are jet-like while BB events are more spherical. The reconstructed mo- 
menta of final state particles is used to form various shape variables (e.g. thrust angle, 
Fox- Wolfram moments, etc.) in order to categorize each event. We follow the scheme de- 
fined in Ref. [T3| that combines seven event shape variables into a Fisher discriminant to 
suppress continuum background. 

Probability density functions (PDFs) for the Fisher discriminant and the cosine of the 
angle between the B flight direction and the beam direction in the T(4S) rest frame are 
combined to form the signal (background) likelihood C s n,) ■ The signal PDFs are determined 
from signal MC simulation; the background PDFs are obtained from the side-band data with 
Mb c < 5.26 GeV/c 2 . We require the likelihood ratio 1Z = C S /(C S + Cb) to be greater than 
0.4 for both ppK + and AAK + modes. These selection criteria suppress approximately 69% 
(66%) of the background while retaining 92% (91%) of the signal for the ppK + (AAK + ) 
mode. In this study only one B candidate is allowed per event. If there are multiple 
B candidates in one event, we select the one with the best x 2 value from the vertex fit. 
Multiple B candidates are found in less than 2% (5%) of events for the ppK + (AAK + ) 
mode. 




M - (GeV/c 2 ) 



FIG. 1: B yield versus M p p. The inset shows the J /ip mass region. (green dots, red dots, and solid 
line represent fitting background shape, signal shape, and combined result, respectively) 
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We use an unbinned likelihood fit to estimate the B yield: 

e -(N s +N b ) N 

L = YliN^M^AE,) + N b P b (M hCi , A^)], 

^ v • i=i 

where P s (Pb) denotes the signal (background) PDF, N is the number of events in the fit, 
and N s and N b are fit parameters. For the signal PDF, we use a Gaussian in Mb c and a 
double Gaussian in AE. We fix the parameters of these functions to values determined 
by MC simulation [18]. Background shapes are fixed from fitting to sideband events in 
the region: 3.14 GeV/c 2 < M pp < 3.34 GeV/c 2 . The M^ c background is modelled using a 
parametrization first used by the ARGUS collaboration, f(M^ c ) oc M^ c \/1 — x 2 exp[— £(1 — 
x 2 )], where x is defined as M hc /E hearn and £ is a fixed value. The AE background shape is 
modeled by a first order polynomial. 

As the mass resolution of M pp (M AA ) is about 10 MeV/c 2 , we determine the B yield as 
a function of M pp (M AA ) from 1.85 GeV/c 2 to 4.5 GeV/c 2 in 10 MeV/c 2 bins. The result is 
shown in Fig.^(Fig.EJ). There are clear 7] c and J/ip peaks in the mass spectrum. A fit to the 
data is shown in the inset. We use a Breit-Wigner function for the r] c peak, a Gaussian for 
the J/ip peak and a line for the non-charmonium background. The background is negligible. 
We define the r/ c signal region as 2.94 GeV/c 2 < M pP < 3.02 GeV/c 2 and the J/ip signal 
region as 3.06 GeV/c 2 < M pp < 3.14 GeV/c 2 . The measured B yield is 329 ± 19 (195 ± 15) 
for B + — > J/ipK + , J/ip — > pp (B + — > r] c K + , r\ c — > pp). We use a phase space MC sample to 
determine the efficiency in this M pp range. The obtained efficiency is 38% (36%) for B + — > 
J/ipK + , J/ip — > pp (B + — > i] c K + , i] c — > pp). The measured branching fractions for charmonia 
decaying into pp are B(r) c -> pp) = (1.58 ± 0.12(stat) ± 0.22(syst) ± 0.47(PDG)) x 10~ 3 and 
B(J/i/> -> pp) = (2.24 ± 0.13{stat) ± 0.31{syst) ± O.Ol(PDG)) x 10" 3 , where the last errors 
are related to the uncertainty in the current world average values of the branching ratios 
B(B + -> r] c K + ) and B(B + -> J/^jK+). 

We study the proton angular distribution in the helicity frame of the J/ip. 9x is defined 
as the angle between the proton flight direction and the direction opposite to the flight of 
the kaon in the J/ip rest frame. The angular distribution is parameterized as P(cos$x) = 
(1 + «b cos 2 6 x )/{2 + 2/3as). Note that a# determined from longitudinally polarized J/ip 
is related to the a determined from transversely polarized J/ip by a# = r^k [111- From 
previous measurements of a, the expectation of is — 0.80±0.04. We modify the likelihood 
function to 

e -(Ns+N b ) N 

L = Y[[N s P s (M hCi ,AE l )e(cos6 x )P(cos6 x ) + N b P b (M hCi ,AE l ,cos6 x )], 

where e(cos#x) is the normalized efficiency function. The observed distribution of e(cos#x) 
is flat. We assume there is no correlation between Mb c , AE and 9 X based on a study with 
MC data. The background PDF, including cos#x> is determined from M pp sideband data. 
The distribution of cos^x for J/ip candidates and fit result in whole M hc ,AE region is 
shown in Fig. El cxb is determined to be —0.54 ± 0.14. 

As a cross check, we fit the 1 + a^cos 2 ^ parametrization to the efficiency corrected 
B yield as a function of cos#x- The results of the fit are shown in Fig. |H1 The value of 
as obtained from the fit is —0.46 ± 0.16 with x 2 /d.o.f. = 1.1. This is consistent with the 
likelihood result from a toy MC study, where samples with the same number of events as 
our data are generated to check the as difference between the likelihood and \ 2 methods. 
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FIG. 2: Likelihood fit of the J/ip — > pp helicity angle distribution. The blue solid, red solid, and 
dashed line represent the fit results, the signal shape, and the background shape, respectively. 



The otB difference follows a Gaussian distribution with a width of approximately 0.05. We 
also apply both the likelihood and \ 2 methods to a J/ip — > ji + fi~ event sample. The result 
is shown as an inset in Fig. El It is in excellent agreement with theoretical prediction which 
has a sin 2 9 X shape (a B = -0.999 ± 0.003). 

We study the systematic error of «g by varying the value of various selection cuts and 
parameters of PDFs to check for trends in the value of a B - This relation is smooth and can 
be fitted to a line. We then quote the change ina B along the line between the selected point 
and the far end of the tested region as a systematic error. Note that this is a conservative 
estimation, since the statistical fluctuation of this data set also contributes to changes in 
a.B- We assign a systematic error of 0.08 for the 1Z selection, 0.06 for PID selection, and 0.02 
for fitting PDFs. Other systematic errors are negligible. To be conservative, we also quote 
the observed difference between the likelihood method and the x 2 method as a systematic 
error. The total systematic uncertainty in a B is 0.13. 

The baryon-antibaryon mass spectrum from B + — > AAK + decays is shown in Fig. 0] 
Similar structures are seen in the mass spectrum of B + — > AAf^ + decays as were seen in 
the B + — > ppK + mass spectrum. There are several complicating factors in the analysis of 
B + — > hAK + decays, relative to B + — > ppK + decays. The slow pion from A decays has 
a low detection efficiency. This causes the A reconstruction efficiency to be non-uniform in 
the polar angle (6 P ) of the secondary decay proton in the A helicity frame, and is correlated 
with cos#x, where X refers to the A. The likelihood function becomes 

£ -(N s +N b ) N 

L = — Y[[N s P s (M hCt , AEi)e(cos 6 X , cos 8 P , cos 6 p )P(cos6 x , cos 6 P , cos 8 P ) 

i=l 

+N b P b (M hCt , AEi, cos9 x , cos 6^, cos 6^)], 
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FIG. 3: J ftp — > pp helicity angle distribution. The dashed line shows the x 2 fit result for B events 
of B + — > J/ipK + , J/ip — > pp. The inset shows the x 2 fit result for S yield of S + — > J/ipK + , J/ip — > 



where e(cos#x, cos# p , cos#p) is determined by a huge signal MC sample with 4 x 10 6 events. 
The background PDF is determined from M AA sideband data in the region 3.14 GeV/c 2 
< M AA < 3.54 GeV/c 2 . The value of a B obtained from the fit is -0.63 ± 0.46 ± 0.27, where 
the systematic error is determined from the same procedure used for J/ip —> pp decays. 

The M bc distribution (with \AE\ < 0.05 GeV) and the AE distribution (with M bc > 5.27 
GeV/c 2 ) for B + — > i] c K + ,i] c — ► AA decays are shown in Fig. |SJ rj c signal peaks are visible in 
the Mbc and AE distributions. The yield from the fit is 19. 5~^'\ with a statistical significance 

of 8.1 standard deviations. The significance is defined as \J — 21n(L /L max ), where L and 
L max are the likelihood values returned by the fit with signal yield fixed at zero and its best fit 
value, respectively. The fit yield is consistent with that obtained in Fig. HJ We estimate the 
branching fraction from the ratio of the efficiency corrected yield of r] c —>■ AA and i] c — > pp. 
The result is B(r] c -»■ AA) = (0.87l^f (stat) ± O.U(syst) ± 0.27(PDG)) x 10~ 3 , where the 
last error is associated with the world average value for B(r] c — > pp) . We apply the same 
procedure to obtain B{J/i/> -> AA) = (2.00±° H(stat) ± 0M(syst) ± 0.08(PDG)) x 10" 3 . 

Systematic uncertainties are studied using high statistics control samples. For proton 
identification, we use aA^ pn~ sample, while for K/n identification we use a D* + — ► D°7T + , 
D° — > K~7i + sample. Tracking efficiency is studied with fully and partially reconstructed D* 
samples. The 1Z continuum suppression uncertainty is studied with b — >• c control samples 
with similar final states. For A reconstruction, we have an additional uncertainty on the 
efficiency for detecting tracks away from the IP. The size of this uncertainty is determined 
from the difference between A proper time distributions in data and MC simulation. Based 
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FIG. 4: B yield versus M AA . The inset shows the rj c mass region. (green dots, red dots, and solid 
line represent fitting background shape, signal shape, and combined result, respectively) 
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FIG. 5: AE and M^ c distribution of B + — » rj c K + ,rj c — > AA candidates. The blue solid, red solid, 
and dashed line represent the fit results, the signal shape, and the background shape, respectively. 



on these studies, we assign a 1% error for each track, 2% for each proton identification, 1% 
for each kaon/pion identification, an additional 3% for A reconstruction and 3% for the TZ 
selection. 

The systematic uncertainty in the fit yield is studied by varying the parameters of the 
signal and background PDFs and is approximately 5%. The MC statistical uncertainty and 
modeling contributes a 5% error. The error on the number of BB pairs is determined to be 
1%, where the branching fractions of T(4S) to neutral and charged BB pairs are assumed to 
be equal. Although the background in the M pP and M AA spectra appear negligible, we forced 
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the B yield to be positive and re-fit the spectra. The feed-down background is estimated to 
be 8% and 12% for the pp and AA modes, respectively. 

To produce a combined systematic error, the correlated errors are added linearly and then 
combined with the uncorrelated errors in quadrature. The total systematic uncertainties are 
14% and 17% for the ppK + , and AAK + modes, respectively. 

In summary, using 386 x 10 6 BB events, we measure the branching fractions of J /if) — > pp, 
Vc ~^ PP, J/ip AA and T) c — > AA from B + — > ppK + and B + — > AAK + decays. We report 
the first observation of r] c -> AA decays, with B(r] c -> AA) = (0.87l[]Jf± 0.14 ±0.27) x 10^ 3 . 
We also measure the parameter as for baryonic J/tp decays. The measured values are 
-0.54 ±0.14 ±0.13 and -0.63 ±0.46 ±0.27 for J/tp -> pp and J/^i -> AA, respectively. The 
above measurements are in agreement with the current world average values as shown in 
TABLE [I] The 5-factories will rapidly accumulate charmonia decays in the coming years, 
enabling more accurate cross checks. 



TABLE I: List of a in previous experiments 



Coll./Mode 


J ftp — ► pp 


J /if) -> AA 


Markl 


1.45 ±0.56 




Mark2 


0.61 ±0.23 


0.72 ±0.36 


Mark3 


0.58 ±0.14 




DASP 


1.70 ± 1.70 




DM2 


0.62 ±0.11 


0.62 ±0.22 


BES 


0.68 ± 0.06 


0.52 ±0.35 


world average 


0.66 ± 0.05(a B = -0.80 ± 0.04) 


0.62 ± 0.17(a B = -0.77 ± 0.13) 
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